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The neutrino-dominated accretion flow (NDAF) with accretion rates M - 0.01-10 A^qS ' 
is a plausible candidate for the central engine of gamma-ray bursts (GRBs). This hyperaccre- 
tion disk is optically thin to neutrinos in the radial direction, therefore the neutrinos produced 
at one radius can travel for a long distance in the disk. Those neutrinos can thus be absorbed 
with certain probability by the disk matter at the other radius and heat the disk there. The 
effect of this "global neutrino heating" has been ignored in previous works and is the focus of 
this paper. We find that around the "ignition" radius ^gn, the global neutrino heating rate could 
be comparable to or even larger than the local viscous heating rate thus must be an important 
process. Two possible consequences are in order if the "global neutrino heating" is taken into 
account: i) the temperature of the disk is slightly raised and the "ignition" radius rign slightly 
shifts to a larger radius, both lead to the increasing of the total neutrino flux; ii) what is more 
interesting is that, the temperature of the ADAF just beyond rign may be raised above the 
virial temperature thus the accretion will be suppressed. In this case, the activity of the black 
hole is expected to oscillate between an active and inactive phases. The timescale of the active 
phases is estimated to be ~ 1 second. If the timescale of the inactive phase is comparable to 
or less than this value, this intermittent activity may explain the slow variability component 
of the GRBs. Self-consistent global calculations of NDAFs with the "global neutrino heating" 
included are required in the future to more precisely evaluate this effect. 
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1 INTRODUCTION 

The accretion disks around black hol es with rather high accretion rates M = 0.01 — 10 M^s'' ar e plausible candidates for the central 
engine of gamma-ray bursts (GRBs) i Naravan et al. 1992 : Popham et al. 19991 : Naravan et all 2001). A large amount of ne utrinos can be 
produced in such hyperaccretion flows. They will then collide, annihilate, and produce a fireball i Meszaros & Rees 2000h . which is one 
popular model for GRBs. In addition to the fireball model, another type of models for the central engine of GRBs is magnetically dominated 
ones. Recent Fermi observations se em to suggest that the latte r is more favored jZhang & Pe'ei l2009l) . The leading magnetic model is the 
Blandford-Znajek (BZ) mechanism i Blandford & Znai"e5ll977 ). which requires a large-scale open magnetic field connecting the black hole 
and an external astrophysical load. This model also tends to generate a contin uous jet. It is shown in this case that the energy dissipation via 
magnetized shoc ks in the jet is not eff icient enough to explain observations i Naravan et al.l 201 1 ). The mag netic turbulent reconnection is 
proposed instead ( Zhang & Yanlboil ). Stimulated by the observations of episodic jets. Yuan & Zhang ( 2012 ) recently proposed a magnetic 
episodic jet model for GRBs. Different from the BZ mechanism, in this model, the jet is intrinsically episodic, i.e., in the form of discrete 
magnetized blobs. Reconnection can be easily triggered when they co llide and GRBs are then natu rally produced. 

In this paper, we focus on the fireball model. The calculations bv lChen & Beloborodovl ( 120071) have shown that this hyperaccretion disk 
has an "ignition" radius r^^^ within which radius the neutrino production become efficient and the disk is neutrino-cooling dominated; while 
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outside which radius the accretion flow is advection dominated. It has been found that except in the very inner part of the disk (i.e. r < 10 r,, 
where = IGMjc- is the Schwarzschild radius of the black hole), the disk is optically thin to neutrinos in the radial direction, therefore the 
neutrinos produced at one radius can travel for a long distance in the disk before they are absorbed or scattered. These neutrinos can thus be 
absorbed with certain probability by the disk matter at the other radius and heat the disk there, as we will show below in this paper. We call 
this effect "global neutrino heating". This process has been ignored in previous works. In this paper we focus on this "global neutrino heating" 
in the neutrino cooling accretion flow, we find that around the "ignition" radius r^^, the global neutrino heating rate could be comparable to 
or even stronger than the local viscous heating rate thus is non-negligible. We note th at similar "global Compton s cattering" effect of photons 
has been studied in the context of optically thin advection-dominated accretion flow ( I Yuan. Xie & Ostrikeiil2009n . 

In section 2, we briefly review the accretion disk model that we are considering. Then in section 3 we show how to calculate the global 
neutrino heating rate and the neutrino radiation pressure. Section 4 is devoted to the numerical results. A brief summary and discussion are 
given in section 5. 



2 THE STRUCTURE OF ACCRETION FLOW AND THE PRODUCTION OF NEUTRINOS 



Consider a stellar-size Schwarzschild black hole (throughout this paper we set the black hole mass M = 3 Mg and the spin a = 0) 
surrounded by a hyper- accreting accretion flow with a typical accretion rate M = 0.1 - 1.0 MgS"'. Within the ignition radius r^^, the 
neutrino cooling is important and t he disk can be regarded as a neutrino-dominated accretion flow (NDAF). The disk properties in this region 
can b e approximately described by jPopham et al .Ii999|JGu et al.l2006lJchen & BeloborodovllOOTi : iLiu et al .I2OO7I : iKawanaka & Mineshigd 
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At the outer disk, r > r^^, when the temperature falls down and the mean electron energy becomes lower than (m,, - ■-" ^'-^ 
trino emission then rapidly switches off and the neutrino cooling becomes inefficient. At this region, the disk thickness is reduced and the 
cooling is dominated by advection. The pr operties of the advection-dominated accretion flow (ADAF) can be approximately estimated by 
I Wang &ZhoiJl 19991 ; lYuan & Zhan32012l) 
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Here T, H, p, Q and v are the temperature, half-thickness, density, angular velocity and radial velocity of the accretion flow respectively; a 
is the viscous parameter, M is the mass accretion rate in unit of MgS 



' M is the black hole mass in unit of Mg, and R is the radius r in unit 



of r,. 



The accretion disk is made of a-particles, neutrons, protons, electrons, positrons, photons, neutrinos and antineutrinos. The outer disk 
is mainly constituted of Q--particles, but once T reaches about 10'" K (happened at around lO^r,), the photodisintegration breaks down a- 
particles into neutrons and protons. The parameter X^^^ is introduced to describe the mass fraction of free nucl eons and then 1 — Xnuc is the 
mass fraction of a-particles. The mass fraction X^^^ is found from the equation of nuclear statistical equilibrium i Chen & Beloborodovl 2007 ; 
Meveilll994h 
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where stands for the proton-to-baryon ratio 
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(15) 



The distributions of the electrons and the positrons are described by the Fermi-Dirac distribution 
1 



gEJk^T^,,, + 1 



(16) 



Here rj^ is a dimensionless degeneracy parameter of electrons defined by r}^ = fi^/k^T, where yu^ (-f^^) is the chemical potential of electrons 
(positrons) and = 8.617x10"" MeV K ' is the Boltzmann constant. If is much larger than unity, the electrons are strongly degenerate, 
whereas if rj^ <K 1, the electrons are weakly degenerate and we can ignore the degeneracy. From Eq. (16) we can obtain the number densities 
of electrons and positrons 

n = — ^ 



-/dp , 



Since the disk matter is neutral, the charge neutrality requires 



(17) 



(18) 



Adopting the conclusion draw bv lChen & Beloborodovl ( 120071) that electrons and positrons are neither nondegenerate nor strongly degenerate 
at all radii, we simply choose j/^ = 1 at all radii and calculate 



(19) 



By taking this simplification, we avoid the iteration in our calculations. Note that at the outer disk the matter is dominated by a-particles 
therefore = 0.5 should be taken as a boundaiy condition. In our calculation we simply set = 0.5 as the upper bound. 

It has been found in previous works that for the hyperaccretion NDAF the most important heating and cooling processes are the viscosity, 
advection and neutrino cooling. The vertically integrated viscous heating rate (over a half-thickness H) is given by 



-(M.MQ.r')n^ 
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In the inner NDAF region, t he advective energy transport rate Q^^^ can be approximately given by ( Naravan & Yilll994 lAbramowicz et al 
I995ljDi Matteo et alj20oi) 
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where Z = 2pH is the surface density of the disk, s denotes the specific entropy, ^ is taken to be constant and equal to 1 and a,. = 4cr^/c is 
the radiation constant. The three terms in the brack ets are the entropy de nsity of radiation, neutrinos, and gas, respectively. And in the outer 
ADAF region, 2", can be approximately given by JWang&Zhoulll999l) 
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where ^ = (4 - 0.75/3)7^ - (12 - 10.5/5)7,,, /? is the ratio of gas to total pressure, = dlnp/dln r, and 7^ = din T/ In r. In this paper we take 
yS = 0, since the radiation pressure dominates over the gas; thus f = 47 - 127^ =^ 1.5. 



In the NDAF, the neutrino or antineutrino production is dominated by the pair capture of electron and positron 



& n + — > p + Vp, which is more familiarly named as the URCA process dKohri & MineshigelEooi : 
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200a). In the case that is of the order of several MeV {E^ <sc m^^ c ), the cross sections of the two processes of the URCA process are 



approximately the same 

G|a + 3£i) 



(23) 



where Gp = 2.302 x 10"^- cm ■ MeV"' is the Fermi coupling constant and « 1.27 is the nucleon axial weak charge. Here E^^-^ is the 
neutrino or antineutrino energy in the final state and we have E^ x E^_ - Q, E^, x E^+ + Q with Q = - mp)c^ ~ 1.29 MeV according to 
the energy conservation. Taking into account the distributions in Eq. (16), we are able to write down the corresponding neutrino cooling rate 
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The emitted monochromatic neutrino (antineutrino) luminosity through the URCA process from a shell at r with thickness dr and height 
H{r) are given by 



Gl(l + 3gl) (P + Q)J(E^ + Qy -m^x'^El 



Gl{\+3g\) (£„ - Q) J(E- - 2)2 - mjc*El 
AL,(E,,r) = lnrAmr-)n„ '^^^^^l' .^av^Jt..^^ . (28) 

And the total luminosity of neutrinos (antineutrinos) from the shell dr are given by 

Kir) = d£.. , (29) 
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L,(r) = I ^^d£, . (30) 



In the NDAF region, the energy loss by neutrino cooling dominates over the energy advection, therefore Q*.^ « . While in the ADAF 
region, we have Q^,.^ a 2^^,,, i.e., the viscous heating is balanced by the advection cooling. 



3 GLOBAL NEUTRINO INTERACTION 



When travelling inside the disk, the emitted neutrinos and antineutrinos may interact with the accretion flow matter and exchange energy 
and momentum with them. We firstly introduce the main interaction channels and the correspondin g energy and momentum loss rate here 
I Burrows & ThomDsonl2002 : iKohri & Mineshigell2002lJchen & BeloborodovlboOTlJltoh et al.ll996h . Since the electron neutrinos (produced 
via the URCA process) are dominated over neutrinos of other two flavors, we consider in this section only the interactions of electron 
neutrinos/antineutrinos with the matter. 

• neutrino absoiption by nucleons: + ;j — > + p and + p ^ + n 

This is the dominant interaction process in the NDAF region since there are plenty of nucleons there. If the energies of neutrinos and 
antineutrinos are relatively low (~ MeV) and <c 'k,/'^ is satisfied, the cross sections of above two processes are approximately the same: 



^.„(£.) - ^^^^-^^E%, « 9.85 X 10-^^4^ 



(31) 

where E^^^-^ is in unit of MeV. Note that, the reaction + p « has a low threshold E^, > 1.806 MeV a m^, + hi,, - m^. When neutrinos 

or antineutrinos are captured by the nucleons, all the energy and the momentum carried by neutrinos will be transferred to the flow matter. 
• neutrino-electron/positron elastic scattering: v + e", e+ — » v + e", e* 

This is the dominant interaction process in the ADAF region since there is little nucleons there. The processes of electron neutrinos 
scattering on electrons or positrons are p roceeded via not only the n eutral current interaction but also charged current interaction. The 
corresponding cross sections are given by (iBurrows & ThompsoijEooi) 
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where 6^ is the W einberg angle and s i n 6^ x 0.23. The average energy transfer from neutrino to electron (positron) is found to be 
(^v(v) ~ ^kgT)/2 bv lTubbs & Schramml 1 119751) . We introduce the parameter 9 to denote the average neutrino energy loss rate in neutrino- 
electron/positron scattering: 9 = {E"^^,^ - E°"l^) / E'"^-^ ^ 1/2 - 2k^T / E^^-y Since the typical energy of neutrinos Ey(y) is roughly equal to the 
typical energy of electron in the radius (kT^) where neutrinos originate, this scattering always heats the disk matter, a s in the NDAF region 
menti oned above. This is quite different from the case of global photons scattering in the case of optically thin ADAF JYuan. Xie & Ostrikei 
20091) . In that case, the energy of photons can be larger or smaller than the average energy of electrons thus the scattering can play a heating 
or cooling role. As for the momentum transfer, the mean scattering angle of this process is quite large (60° ~ 90°), which means that even for 
relatively small energy transfer there can be a large momentum transfer in this scattering. In our numerical calculations, we simply choose 
the momentum transfer cross sections to be cr"^ = 2a^,J3. 

• neutrino-baryon elastic scattering: y + p, «, He — > v + p, n. He 

Neutrinos and antineutrinos interact elastically with nucleons («, p) or the u-particles via the neutral current interaction. The corresponding 
cross sections are the same for either neutrinos or antineutrinos and are approximately given by 
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Note that, the MeV neutrinos scatter from nuclei or nucleons almost coherently, which means neutrinos do not loss energy in these scattering. 
However, the mean scattering angles of these processes are large, which means there are also large momentum transfer in the scattering. In 
our calculations, we also simply choose the momentum transfer cross sections to be cr^ „ „ = 2a „ ^^/3 and ignore the energy transfer. 

In order to calculate the "global" neutrino heating rate at a given radius r, we nee d to know the received neutrino spectrum at r. The 
received neutrino (antineutrino) spectrum at r emitted by the flow inside of r is given by jYuan. Xie & Ostrikeill2009l;|Park & OstrikeijboO?!) 



1 dL^,„(£^,,), r') 
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where r.,, = 3r^ is the inner edge of the NDAF, dL^^-^(E^^-y r') has been given in Eqs. (27) and (28) and t^^-^ is the neutrino (antineutrino) 
optical depth from r' to r 



with ((Tn)j,|^-j are functions of r" and r' 

{an\ = tr,„(£;)< + a ^,{E'X, + a/E[X + <^.(KK + 'r„(K)(n';^ + <_) , 
(o-n), = cr,p(£;);i;' + cr i,(E',)n'i, + a j,iE[)n';, + (r„{E',X; + cr-JE'„)(n'^- + «;.'_) , 

where cr^„ ^ „ ^ can be calculated by using Eqs. (31) - (36), in which E'^^-^ is the neutrino/antineutrino energy emitted at the radius r' 
in unit of MeV, T" is the temperature at r" and n"^ ^ ^ are the number densities of cor responding particles at the r adius r". The neutrino 
(antineutrino) spectrum received at r emitted by the flow outside of the rdius r is given bv J Yuan. Xie & Ostrikej2009HPark & Ostrikeij2007l) 
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lO^Tj in our calculations. The results are not sensitive to the exact value of r,„„. The total neutrino (antineutrino) spectrum 



received at r is the sum of F"-^ and F""ly Then the global neutrino heating rate can be calculated by 
q.iobai = f [pyiE,., r) + F^iE,, r)] [aJE,)n„E, + crJE^.)(n^^ + )6E,) dE, 

+ f [f^"(£„ r) + E^"(E„ r)] [o-,^,(E,)n„E, + o-,^{E,)(n^^ + n^-)9E,) dE, , 
where n^^ are the number densities at the radius r. The vertically integrated global neutrino heating rate is given by 



(42) 



H 



(43) 



Q global ~ 'J global ' 

Besides the energy transfer from neutrinos to the flow matter, the emitted neutrinos/antineutrinos also transfer momenta to the flow 
matter when interacting with them, i.e., the radiated neutrinos exert pressure upon the flow matter. The average force suffered by the disk 
matter can be approximately evaluated by 
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where neutrinos are regarded to be massless and n,^, = + + n^. This force is placed against the gravitational force of the central black 
hole 



(M ■ Mg)mp 



(45) 



where G„ is the gravitational constant. 



4 NUMERICAL RESULTS 

We consider in this paper four accretion flow models: (a) M = 0.1 MqS"', o- = 0.1; (b) M = 1.0 MqS"',^ = 0.1; (c) M = 0.1 MgS~',g = 
0.01;and(d)M = 1.0 M^s^^a = 0.5. We calculate the global neutrino heating rate numericallv. Following the results in lChen & Beloborodov 
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Figure 1. The heating and coohng rates Q^,-^ (dashed), Q^i^ij^^i (sohd), (dotted) and 2,, (dot-dashed) as functions of radius r for Model (a) (top-left; 
M = 0.1 MqS"', Q- = 0.1, K 60 r,). Model (b) (top-right; M = 1.0 M^s'' . a = 0.1, r^^ a; 200 r,). Model (c) (bottom-left; M = 0.1 M^s"', a = 0.01, 
''ien ~ Model (d) (bottom-right; M = 1.0 MqS"', a = 0.5, r^^ » 40 r,), respectively. 



1 20071) . we set r^^ = 60, 200, 450 and 40 r, in Models (a), (b), (c), and (d), respectively. The four plots in Fig. 1 show the viscous heating rate 
(2,!,, (dashed), global neutrino heating rate (solid), advection cooling rate 2^^^, (dotted), and neutrino cooling rate Q^, (dot-dashed) for 

Models (a)-(d), respectively. Some quantities are not continuous across r^^, this is because the analytical solutions of Eqs. (I)-(IO) are not 
precise. 

From Fig. 1 we can obtain the following results. 

• At the innermost region of the NDAF, the global neutrino heating rate is always much lower than the local viscous heating rate. 

• However, with the increase of the radius, the ratio 2g;o4o//2w, increases, i.e., the global neutrino heating becomes more and more 
important. Physically, this is because the two terms have different scaling with radius. From Eq. ( 1201 we see that Q'^^.^ oc The radial 
scaling of Qgi^ij^i is not easy to estimate. From Eq. ( 1431 it is roughly proportional to F^^-^n^^ ^r. Here F^^-^ is the received neutrino flux at r, its 
radial scaling must be flatter than because of the accumulation of neutrinos from all the radii within r. The number density ;i„ ^ oc r"-^^ 
from Eq. (|4j. Thus it is quite possible to have the scaling of Q*i„i,^i flatter than r"^. 

• It is interesting to note that around the "ignition" radius r^^^^, the global neutrino heating rate becomes comparable to the viscous 
heating rate Q*.^ for Models (a), (b), and (d). For Model (d) the global heating even becomes larger than the viscous heating in the ADAF 
region. This result is especially of interest to us, as we will discuss below. 
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of the protons in the same accretion flow. 



• When the radius is larger than r^^, the global heating rate rapidly decreases with increasing radius. This is because the temperature 
decreases, thus ;i„ ^ quickly decreases. 

• By comparing the four plots in Fig. 1, we can see that the ratio 2g,„i,„,/2J,j around r^^ is more sensitive to a than to M. A higher a, 
rather than a higher M, can more easily result in a larger ratio. Physically, this is because the produced total neutrino flux in the NDAF is 
fully determined by M. But a larger a implies a smaller density thus smaller radial optical depth, which then results in a stronger neutrino 
flux received at r. 

We can see that in some case such as Model (d), the global neutrino heating is very significant. It is even larger than the viscous heating 
close to r^j^ in the ADAF region. In this paper we only focus on the evaluation of the significance of the global neutrino heating, but do 
not self-consistently calculate the dynamical structure and radiation of the hyperaccretion flow after the global neutrino effect is taken into 
account. However, two consequences can be expected: 1) the temperature of the disk at around r ^ will be slightly raised and the ignition 
radius r^^^ will become larger. In this case the neutrino emission will be enhanced since the efficiency of the UCAR process is very sensitive 
to the temperature; 2) more importantly, if the global heating effect is important enough in the ADAF region, since the temperature of the flow 
in the ADAF is already close to the virial temperature, it is possible that the temperature of the ADAF will be higher than the virial one due 
to this global heating. In this case the accretion rate will be highly suppressed after the matter inside of r-^^^ falls onto the black hole horizon. 
Only when the heated flow cools down or the unheated flow located at large radius comes in, does the accretion rate reco ver. This means tha t 



the black hole activity will oscillate between an active and an inactive phases. Su ch kind of behavio r was first found bv lCowie et al.l i\91 



and later was used to explain the intermittent activity of compact radio sources (lYuan & Lill201 ih . The duration of the inactive phases is 
determined by the cooling timescale of the heated flow, or the accretion timescale of the unheated flow located at large radius. The duration 
of the active phase is roughly determined by the accretion timescale of the NDAF at r^^, 

~ r/v ~ 5 X M'-^a-'-\r-^Jrf^ ~ 1.1(M/3Mq)' ^(a/O.l)-' 2(^_^^^yiOO,-,)°« s (46) 



Such a timescale is roughly equal to the "slow" variability timescale of the GRBs i Gao et aill2012l) . If the duration of the inactive phas 



lase IS 



comparable to or less than Tad, the global neutrino heating mechanism may explain the origin of the "slow" variability component of GRBs; 
otherwise this may imply a challenge to the NDAF model. A more exact global calculation with the global neutrino heating effect included 
is badly needed to obtain the self-consistent solution of the heated flow and thus make a final conclusion. 

In addition to the energetics, the average neutrino radiation pressure on the disk matter is also calculated numerically. The results for 
Model (a) are illustrated in Fig. 2. One can see that the average neutrino radiation pressure Ty (solid) is much lower than the gravitational 
force of the protons (dashed). We conclude that the neutrino radiation pressure in this hyperaccretion flow can be safely neglected. 



5 SUMMARY AND DISCUSSIONS 

In a NDAF model, most of the neutrinos are produced at the inner region of the accretion flow. Since the optical depth for neutrino is small 
in the radial direction, these neutrinos will be able to travel for a long distance and transfer their energy and momentum to the accretion 
matter at large radii. This "global neutrino interaction" effect has been neglected in previous works and its significance is the focus of the 
present paper. Our numerical analysis shows that if a is not too small, o- > 0. 1 , the global neutrino heating rate will be comparable to or even 
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larger than the local viscous heating rate at r > r-^^^^ (here r_,^ is the "ignition" radius, i.e., the boundary between the inner NDAF and the 
outer ADAF). Consequently, the temperature of the ADAF may become higher than the virial one. This implies that the flow outside of r^^ 
can't be accreted and the accretion will be highly suppressed once the accretion flow inside of r^^^ is accreted. We found that the accretion 
timescale of NDAF at r^^ is typically 1 second. This means that the NDAF model will not be able to explain the long GRBs. We would 
like to emphasize, however, that our calculations are based on the simple analytical solutions of NDAF and ADAF and are "local", therefore 
the results should only be taken as preliminary. To get a more precise answer to the significance of the global neutrino interaction effect, a 
self-consistent and global calculation with this effect included is required. The readers can refer to Yuan, Xie & Ostriker (2009) for the details 
of such a calculation. Here we only briefly review the approach. The first step is to get the global solution of the set of difi"erential equations 
describing the superaccretion flows (e.g., Popham et al. 1999). No global neutrino heating term is included in this step. Then we follow a 
similar approach as presented in this paper to calculate the global neutrino heating rate at each radius ^'j/oia/C'')- Then we include this term 
in the energy equations of the accretion flow and try to solve their global solution again. Based on this new global solution we can calculate 
^giobai^'^ again and compare with the result obtained in the last step. If they are not equal to each other, we replace 9j,„,„,('') with the new one 
and repeat the above procedure until convergence is achieved. 

In addition to the "local" nature of our calculation mentioned above, another caveat of our work is that we do not consider outflow 
in both NDAFs and ADAFs. The rad iation-megnetohydrodynamical numerical simulations to optically thick ADAFs show the existence of 
outflow jOhsuga & MineshigeluOl ih . The mechanism is that the radiation force is found to exceed the gravitational force at the surface of 
the flow. We are not aware of any theoretical work focused on the st udy of possible outflow in NDAFs. However, we note that the "micro 
Blandford-Payne" mechanism recently proposed by I Yuan et al.l ( 120 12h to explain the origin of outflow in an optically thin ADAF should work 
for any accretion flow in which the magnetic turbulence associated with the magnetorotational instability is the mechanism of transporting 
angular momentum. This mechanism is veiy similar to the Blandford & Payne one in the sense that the magnetic centrifugal force drives 
the outflow. The difference is that no large-scale open field is required here (see also Bai & Stone (2013, preprint) for the shearing box 
simulation analysis). Th e initial sim ulations to the radiatively efficient standard thin disk, which is physically similar to the NDAF, do show 
the existence of outflow jOhsuga & Mineshige 201 li). However, all these are our speculation since the parameters of accretion flows here are 
largely different from the works mentioned above. It is hard to estimate the effect of outflow on the significance of global neutrino interaction. 
On one hand, the inward decrease of accretion rate will make the neutrino production in the inner region of NDAF weaker in relative to the 
viscous heating rate at r-^^^. This will make the global neutrino effect less important. But on the other hand, the inward decrease of accretion 
rate means that the radial optical depth will become smaller which enhances the received neutrino flux at r-^^^. The final result requires detailed 
numerical calculation. 
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